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Deducing the asymptotic normalization constant 
of the 2 + subthreshold state in 16 from 12 C + a elastic scattering 
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R- matrix analyses of the 12 C + a elastic-scattering phase shifts deduced from a recent 
high-precision measurement of the differential cross sections are performed. The i = 
phase shifts constrain the R-matrix radius a around 5.85 fm, while the £ = 2 phase shifts 
lead to a strong constrain neither on a nor on the asymptotic normalization constant C of 
the 2 + subthreshold state (except for a loose upper limit). This contradicts previous R- 
matrix analyses of the 12 C + a elastic scattering and explains the incompatibility between 
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values of C obtained in these analyses. 



1. Introduction 



The 12 C (a, 7) 16 O capture reaction is very important in nuclear astrophysics: combined 
with the triple-a reaction, it determines the rate of carbon and oxygen resulting from the 
helium-burning phase in red giants. Measuring this reaction rate is a difficult task: due 
to the Coulomb barrier, the cross section is very small (10~ 17 barns at E cia = 300 keV, 
the astrophysically-relevant center-of-mass energy). Hence, a theoretical extrapolation 
from measured higher energies is necessary, which is generally performed with the R- 
matrix formalism. However, this extrapolation is made very unstable by the presence 
of two bound states in the 16 spectrum just below the 12 C + a threshold: a 1~ state 
strongly influences the dominant El capture, while a 2 + state strongly affects the smaller 
E2 component. Despite some controversy [[Q, the El component is generally believed to 
be well constrained by the measurement of the /5-delayed a spectrum of 16 N leaving 
the influence of the 2 + subthreshold state as the dominant source of uncertainty on the 
total reaction rate. 

In the R-matrix formalism, the key parameter is the bound-state asymptotic normal- 
ization constant (ANC), which is related to the reduced width. Since this quantity also 
appears in the description of elastic scattering, it has been proposed to constrain it by 
precise measurement of elastic-scattering differential cross sections. An R-matrix anal- 
ysis of the 2 + phase shifts deduced from an older experiment [ Ej leads to the value 
C = 402 x 10 3 fm" 1 / 2 , but no error bar is provided [ Hj . A recent high-precision mea- 
surement of N = 402 differential cross sections at center-of-mass energies between 2 and 
6.5 MeV, motivated by the astrophysical context described above, has been performed at 
Notre-Dame [Ej. Though the experimental results are compatible with the older experi- 
ment, their R-matrix analysis leads to a totally different ANC value: C = 154 ± 18 x 10 3 
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f m -i/ 2 This discrepancy suggests a problem in these R-matrix analyses, which motivates 
the present work. 

In the next section, I perform new R-matrix analyses of the phase shifts deduced from 
the Notre-Dame experiment and solve the above contradiction. A brief comparison 
with the potential model is then made in the last section. 

2. R-matrix analyses of the 12 C + a elastic-scattering phase shifts 

In the R-matrix description of elastic scattering, a partial-wave decomposition of the dif- 
ferential cross section is performed. The phase shifts of a given partial wave are expressed 
as the sum of a hard-sphere background term with radius a and an R-matrix or "reso- 
nant" term which is parameterized in terms of energies and reduced widths of the states 
in the given partial wave. These states may be actual resonances at positive energies, 
which lead to rapid variations in the elastic-scattering excitation functions and are exper- 
imentally visible, but they may also lie outside the measured energy region and only lead 
to a smooth background influence in the measured region. These background terms may 
correspond either to bound states (negative energy), which have a clear physical meaning, 
or to high-(positive-)energy states, which are generally introduced for phenomenological 
reasons and have a less clear physical meaning. 
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Figure 1. i = phase shifts as a function of the center-of-mass energy: experimental 
points [ |H] , R-matrix fit with a background consisting of (i) a hard-sphere term only 
(dashed line) and (ii) a hard-sphere and a high-energy-state terms (full line). 

Figure ^ shows a fit of the £ = 12 C + a elastic-scattering phase shifts. There is only 
one narrow resonance at 4.9 MeV, superimposed on a smooth background. Let us examine 
the different contributions to that background by first considering a simple hard-sphere 
term. The dashed curve is the best fit, obtained for a = 3.66 fm; it corresponds to a deep 
minimum of \ 2 as a function of a, as shown on Fig. |2l but it is clearly not satisfactory. 
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Figure 2. yX 2 /N on the iV = 402 £ = phase shifts from Ref. [ |Sj for an R-matrix 
fit with a background consisting of (i) a hard-sphere term only (dashed line) and (ii) a 
hard-sphere and a high-energy-state terms (full line), as a function of the R-matrix radius. 

The first way to improve this fit is to consider bound states. I have checked that the two 
+ bound states of 16 O, when added into the fit, do not improve the situation: they are 
too far away from the measured region. In Ref. [Ej, the state is however taken into 
account into the fit and a very small reduced width (and hence a very small ANC) is 
obtained. I believe this small ANC has no physical meaning: it only indicates that the 
fit is not less good when the bound state is not included. The second way to improve the 
fit is to consider a high-energy phenomenological state. The full curves in Figs. ^ and El 
correspond to that case. The fit is now satisfactory; it corresponds to a clear minimum 
of x 2 f° r a — 5.85 fm, which is close to the value of 5.5 fm obtained in Ref. [E]. 

The present analysis thus confirms that the R-matrix radius can be constrained by 
elastic scattering, as claimed in Ref. [E|. The origin of this constrain is mostly the I = 
partial wave: I have checked that the minimum of \ 2 is much less marked for other partial 
waves (see below the case of I = 2). This can be explained by the fact that the background 
phase shifts decrease for increasing angular momentum. The present results also show 
that the value of the radius strongly depends on the presence or absence of a high-energy 
background state. While this radius is physically understandable in the absence of high- 
energy state (3.66 fm is close to the sum of the matter radii of the 12 C and a nuclei), the 
fit is unsatisfactory; on the other hand, with a high-energy state, the fit becomes good 
but the physical meaning of the radius, 5.85 fm, is less clear. 

Let us now go to higher angular momenta: for each partial wave, a high-energy state 
is necessary to get a satisfactory fit of the phase shifts, as for I = 0. The background 
contains thus at least two terms. To extract a bound-state ANC from the phase shifts, one 
has thus to disentangle the bound-state background contribution from a total background 
which consists of three terms (see Fig. 2 of Ref. [Ej). I have shown above that this was 
not possible for i = 0. The same conclusion holds for the 3~ bound state, in contrast 
to what is done in Ref. [|5]: the fit of the phase shifts is not better with the bound 



4 



state. In both cases, this can be understood by the rather low energy of the bound states; 
for subthreshold states, one could expect the situation to be better. For the 1~ state, 
however, no evidence of better fit is obtained in Ref . [ |3] . The general tendency thus 
seems to be that bound states do not make the fit of elastic scattering better. With this 
in mind, the results of Refs. [ |1] and [ |3] for £ = 2 are rather surprising. 




Figure 3. Same as Fig. [T] but for £ = 2. 



I have tried to reproduce these results but did not find any relevant x 2 minimum which 
could lead to a constrain on the ANC. The situation for £ = 2 is actually very similar to 
that for £ = 0, as shown by Fig. El a good fit requires a high-energy background state but 
does not require a bound-state term. The situation is best illustrated on Fig. EJ where the 
X 2 on the phase shifts is plotted as a function of a and C. There is a valley of satisfactory 
values for these parameters, which ranges from C = (no bound state in the fit) to at 
least C = 200 x 10 3 fni^ 1 / 2 . These fits are satisfactory in the sense that they give a good 
graphical fit of the phase shifts (see Fig. |3J) but they are not totally satisfactory from 
the point of view of x 2 theory, which provides me from rigorously defining the acceptable 
values of C. This is probably due to an underestimation of systematic errors on the phase 
shifts of Ref. [ |U] . Figure 0] however clearly shows that the R-matrix fit only provides an 
upper limit on the ANC, strongly depending on a. The conclusion is thus the same for 
all partial waves: an R-matrix analysis of a single partial wave is not able to constrain 
a bound-state ANC, even when the bound state is close to threshold. This shows that 
the disagreeing values of C obtained in Refs. [0] and [E] are unreliable. I think that the 
X 2 minimum obtained in Ref. [ E| is due to too restrictive hypotheses on the R-matrix 
radius and the high-energy state, and that the error bar on C obtained in Ref. neglects 
systematic errors on the data and is therefore much too small. 
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Figure 4. \jx 2 /N on the N = 402 t = 2 phase shifts from Ref. [ [6J for an R-matrix 
fit with a background consisting a hard-sphere, a high-energy-state and a bound-state 
terms, as a function of the R-matrix radius a and the bound-state ANC C. The best fits 
correspond to Jx 2 /N ~ 3.7. 



3. Comparison with the potential model 

The above result is not very surprising from the point of view of the potential model: 
it is well known that, for a given partial wave, different potentials may have bound states 
at the same energy but different ANC, for the same scattering phase shifts at all energies. 
Such phase-equivalent potentials are called Bargmann potentials [IZIIHI and are well known 
in the context of the inverse problem [Ej. Such potentials have been constructed for 12 C + 
a in Ref. [QUI- In the same reference, it is shown that the ANC may however be extracted 
from the elastic-scattering phase shifts, in the framework of the potential model, provided 
several partial waves are considered simultaneously. For 12 C + a, this can be achieved 
by fitting the properties of the 0^ rotational band, to which the 2 + subthreshold state 
belongs. This leads to the value C = 144.5 ± 8.5 x 10 3 fm -1 / 2 . 
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